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Magnetic Order Beyond RKKY in the Classical Kondo Lattice 
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Abstract. - We study the Kondo lattice model of band electrons coupled to classical spins, 
in three dimensions, using a combination of variational calculation and Monte Carlo. We use 
the weak coupling 'RKKY' window and the strong coupling regime as benchmarks, but focus on 
the physically relevant intermediate coupling regime. Even for modest electron-spin coupling the 
phase boundaries move away from the RKKY results, the non interacting Fermi surface no longer 
dictates magnetic order, and weak coupling 'spiral' phases give way to collinear order. We use 
these results to revisit the classic problem of 4/ magnetism and demonstrate how both electronic 
structure and coupling effects beyond RKKY control the magnetism in these materials. 



The Kondo lattice model describes local moments on a 
Ijattice coupled to an electron band. Such local moments 
arise from electron correlation and Hunds coupling in the d 
shells of transition metals or the / shells of rare earths. Al- 
though historically the 'Kondo lattice' arose as the lattice 
version [1] of the Kondo impurity problem, and refers to 
S = 1/2 moments coupled to conduction electrons, there 
are also systems with local electron-spin coupling where 
the moment is due to a spin with 25* ^ 1. In that case 
tTie quantum fluctuations of the local moment, and the 
Kondo effect itself, are not relevant. Such a system can 
be described by a classical Kondo lattice model (CKLM). 
This limit is relevant for a wide variety of materials, e.g, 
t|he manganites [2], where S* = 3/2 moments couple to 
itinerant electrons via Hunds coupling, or 4/ metals [3-6], 
e.g, Gd with S = 7/2, or the Mn based dilute magnetic 
semiconductors [7] where S = 5/2. In some of these mate- 
rials, notably the manganites and the magnetic semicon- 
ductors, the coupling scale is known to be large, while in 
the / metals they have been traditionally treated as being 
weak. 

The CKLM involves the ordering of 'classical' spins, but 
the effective interaction between spins is mediated by elec- 
tron delocalisation and cannot be described by a short 
range model. In fact the major theoretical difficulty in 
analysing these systems is the absence of any simple clas- 
sical spin model. Nevertheless, there are two limits where 
the CKLM is well understood, (a) . When the electron-spin 
coupling is small, one can perturbatively 'integrate out' 
the electrons and obtain the celebrated Ruderman-Kittel- 



Kasuya-Yosida (RKKY) model [8] . The effective spin-spin 
interaction in this limit is oscillatory and long range, con- 
trolled by the free electron susceptibility, Xolq), and the 
magnetic ground state is generally a spiral, (b). When 
the electron-spin coupling is very large compared to the 
kinetic energy, the 'double exchange' (DE) limit, the elec- 
tron spin is 'slaved' to the orientation of the core spin 
and the electronic energy is minimised by a ferromagnetic 
(FM) background [9]. This leads intuitively to a spin po- 
larised ground state. 

In many materials the ratio of coupling to hopping scale 
is > 1, but not quite in the double exchange limit. In that 
case one has to solve the coupled spin-fermion model from 
first principles. Doing so, particularly in three dimensions 
and at finite temperature, has been a challenge. We study 
this problem using a combination of variational calculation 
and full spin-fermion Monte Carlo. 

Our principal results are the following: (i) We are able 
to map out the magnetic ground state all the way from the 
RKKY limit to double exchange, revealing the intricate 
evolution with coupling strength, (ii) We demonstrate 
that the phase boundaries depend sensitively on electronic 
hopping parameters. This is not surprising in the RKKY 
regime, but the dependence at stronger coupling is un- 
known, (iii) We use our results to revisit the classic 4/ 
magnets, widely modelled as RKKY systems, and suggest 
that with increasing 4/ moment, the effective coupling in 
these systems pushes them beyond the RKKY regime. We 
work out the signatures of this 'physics beyond RKKY'. 
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Model: The Koiido lattice model is given by 

H =-Y^ Ujc\^Cj„ ~ l^'^rii - ai.Si (1) 

We will use f = 1 as the nearest neighbour hopping ampli- 
tude, and explore a range of i', the next neighbour hop- 
ping, on a cubic lattice. Changing t' will allow us to ex- 
plore changes in the (bare) Fermi surface, and particle-hole 
asymmetry, is the chemical potential, and J > is the 
local electron-spin coupling. We assume the Si to be clas- 
sical unit vectors, and absorb the magnitude of the core 
spin into J wherever necessary, ffi is the electron spin op- 
erator. Wc work with /i. rather than electron density (n), 
as the control variable so that regimes of phase separation 
(PS) can be detected, and study the magnetic properties 
for varying n, i'/t, J/i, and temperature T /t. 

Although there have been many studies in the 'double 
exchange' {J/t oo) limit [10], the attempts to explore 
the full n^J—T phase diagram have been limited, (a) An 
effective action obtained from the CKLM via gradient ex- 
pansion [11] has been analysed. This mapped out some 
of the commensurate and spiral phases in two dimensions, 
where the phases are fewer. It did not explore the finite 
temperature physics, e.g, the Tc scales, and seems to be 
inaccurate when handling commensurability effects near 
n = 1. (b) The model has been studied within dynamical 
mean field theory [12] (DMFT), and the broad regimes of 
fcrromagnetism, antiferromagnetism (AFM), and incom- 
mensurate order have been mapped out. Unfortunately 
the effective 'single site' character of DMFT does not al- 
low a characterisation of the incommensurate phases and 
misses out on the richness of the phase diagram. The loss 
of information about spatial fluctuations also means that 
critical properties, either in magnetism or transport, can- 
not be correctly captured, (c) An 'equation of motion' 
approach [13] has been employed to study general finite 
S spins coupled to fermions, and results have been ob- 
tained in the classical limit as well. However, except the 
ferro and antiferromagnetic phases other magnetic states 
do not seem to have been explored, {d) The full spin- 
fermion Monte Carlo, using exact diagnolisation, has been 
employed [14] in one and two dimensions but severe size 
limitations prevent access to non trivial ordered states. 

Method: The problem is technically difficult because it 
involves coupled quantum and classical degrees of free- 
dom, and there is in general no equivalent classical spin 
Hamiltonian. The probability distribution for spin con- 
figurations is given by P{S} oc Tr^, j,te~^^ so the 'effec- 
tive Hamiltonian' is i?e//{S} = —^log Tr^ ,^-te~^^ , the 
fermion free energy in an arbitrary background {S^}. It 
cannot be analytically calculated except when J/t <^ 1. 

When J/t ^ 1, the (free) energy calculated perturba- 
tively to 0{J'^) leads to the RKKY spin Hamiltonian [8], 

Hrkky = Ey Ji3^t-^3, where J„- ~ J^xlj and x°j is the 
non local susceptibility of the free (J = 0) electron sys- 
tem. Xii is long range and oscillatory. We will analyse this 



model to understand the weak coupling phases. At strong 
coupling, J/t —^ oo, there is no exact analytic Heff but 
we can construct approximate self consistent models [15] 
of the form H'^^ — — Dijy^l + S,;.Sj, with the Dij 
related to the electronic kinetic energy. Unfortunately, 
when J ~ 0{t) neither the RKKY model nor the DE ap- 
proximation are valid. This regime requires new tools and 
we will use a combination of (i) variational calculation 
(VC) [16] for the magnetic ground state, and (ii) spin- 
fermion Monte Carlo using a 'travelling cluster' approxi- 
mation [17] (TCA-MC) at finite temperature. 

For the variational calculation we choose a simple 
parametrisation [18] for the spin configuration: Siz = a, 
Six = V^ — o? cosq.Ti and Siy = \/\ — o? sinq.ri. This 
encompasses the standard ferromagnet and antiferromag- 
net. as well as planar spiral phases, canted fcrromagncts, 
and A and C type antiferromagnets. For a fixed /x and J 
we compute the electronic energy f (a, q, //) and minimise 
it with respect to a and q. The electronic density at the 
chosen /i is computed on the minimised state. Since the 
magnetic background only mixes electronic states |k, t) 
and |k — q, I) the electronic eigenvalues e''=(k, q) are sim- 
ple, and only an elementary numerical sum is required to 
calculate £{a, q, fi) = J2k ± ^"^O^i q)^(M ^ e='=(k, q)). 

While the VC provides a feeling for the possible ground 
states, it has the limitation that (i) it samples only one 
family of (periodic) functions in arriving at the ground 
state, and (ii) finite temperature properties, e.g, the mag- 
netisation and the critical temperature are not accessible. 
For this we 'anneal' the system towards the equilibrium 
distribution P{S} cx Tr^^cte"'^^ using the TCA based 
Monte Carlo. In this method the acceptance of a spin up- 
date is determined by diagonalising a cluster Hamiltonian 
constructed around the update site, and avoids iterative 
diagonalisation [14] of the full system. We can access sys- 
tem size ^ 10^ using a moving cluster of size 4'^. 

The TCA captures phases with commensurate wavevec- 
tor Q quite accurately, but access to the weak coupling 
incommensurate phases is poor. To get an impression 
of the ordering temperature for these phases we compute 
the energy difference A£{n,J) = j^{£disord - £ord), be- 
tween the ordered state and a fully spin disordered state 
in a large system. Eord is calculated from the variational 
ground state, and Edisord by diagonalising the electron sys- 
tem in a fully spin disordered background on a large lat- 
tice. A£{n, J) is the 'condensation energy' of the ordered 
state, and provides a crude measure of the effective ex- 
change and Tc- Where wc could compare the trend to MC 
data, the agreement was reasonable. 

Results: The results of the variational calculation in the 
'symmetric' {t' = 0) case are shown in Fig.l. We employed 
a grid with upto 40^ k points, and have checked stability 
with respect to grid size. Let us analyse the weak and 
strong coupling regimes first before getting to the more 
complex intermediate coupling regime. 

(i) RKKY limit: The key features for J/t are: 
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Fig. 1: Colour online: Magnetic ground state for the particle- 
hole symmetric model {t' = 0) for varying electron density (n) 
and electron-spin coupling (J). The phases are characterised 
by their ordering wavevector Q, indicated by the colour code 
in the legend to the right, and their net magnetisation a (if 
any). Among the 'commensurate' phases, Q — {0,0,0} is the 
usual ferromagnet, {0, tt, n} and {tt, tt, tt} are antiferromagnets 
with no net magnetisation, while the {0, 0, tt} antiferromagnet 
has a — for J — > but picks up finite magnetisation with 
increasing J. At n = 1 the system is always a Q = {tt, n, n} 
antiferromagnet. The incommensurate phases have ordering 
wavevectors {Qx,Qy,Qz} of which at least one component is 
neither nor tt. For such phases the exact wavevector depends 
on the value of ?i and J. For example, for J ^ the (blue) 
checkerboard region in the left hand corner, to the right of Q = 
{0, 0, tt}, has wavevector Q = {0, Qy,n}, where Qy varies from 
to TT as one moves left to right. The (green) shaded regions in 
the phase diagram, not indicated in the legend, are windows of 
phase separation. No homogeneous phases are allowed in these 
regions. The results in this figure are based on a variational 
calculation using a 20"' k point grid, and cross-checked with 
data on 40^ . 
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Fig. 2: Colour online: The finite temperature phase diagram in 
the particle-hole symmetric case, for various J. Panels (a)-(d) 
show the different ordered phases and their estimated transi- 
tion temperature as we move from the weak coupling to the 
double exchange limit. The legend for the phases is shown 
on the right. The transition temperatures are based either on 
Monte Carlo results (shown as symbols), or the A£ estimate 
(firm lines) described in the text. Notice that the Tc for the 
ferromagnetic, Q — {0,0,0}, phase increases (and saturates) 
with increasing J. At n = 1 the order is at Q = {tt, tt, tt} and 
the corresponding Tc initially increases with increasing J and 
then decreases. Except for Q = {0, 0, 0} and {tt, tt, tt} other 
phases vanish by the time J/t = 10. The Monte Carlo esti- 
mate of ferromagnetic Tc axe shown as circles, while that of 
the antiferromagnet is marked on the n — 1 axis by a square 
symbol. As in Fig.l the (green) shaded regions indicate phase 
separation. 



(i) the occurence of 'commensurate' planar spiral phases, 
with wavenumber Q virhich is {0,0,0}, or {0,0, tt}, etc, 
over finite density windows, (ii) the presence of planar 
spirals with incommensurate Q over certain density in- 
tervals, (iii) the absence of any phase separation, i.e. 
only second order phase boundaries, and (iv) the pres- 
ence of a 'G type', Q = {7r,7r,7r}, antiferromagnet at 
n — 1. Although the magnetic state is obtained from 
the variational calculation, much insight can be gained by 
analysing the H'^^^^y- Since the spin-spin interaction is 
long range it is useful to study the Fourier transformed ver- 
sion H^fjfj^y ^ Eq^qlSqP, where Jq = E,-, ^..e''' ''- 
and Sq = EiS^e^i-^'. The coupling Jq = ,Pxo{(l,n) 
is controlled by the spin susceptibility, xo{^,n), of the 
J = tight binding electron system. For our choice of 
variational state the minimum of H'^-^-^ corresponds to the 
wavevector at which xo (Qj has a maximum. We inde- 
pendently computed Xo(qi'^) a-nd confirmed [19] that the 
wavevector Q{n) obtained from the VC closely matches 



the wavevector q_max{n) of the peak in XoCQi'^)- The ab- 
solute maximum in Xoi^^n) remains at q = {0,0,0}, as 
the electron density is increased from n = 0, and at a crit- 
ical density q^max shifts to {0, 0, tt}. With further increase 
in density qmax evolves through {0, q, tt} to the C type 
{0, TT, tt}, then {q, tt, tt}, and finally the G type AFM with 
{tt, 7r,7r}, where the Fermi surface is nested. The absence 
of 'conical' phases, with finite (a) and a spiral wavevector, 
is consistent with what is known in the RKKY problem. 

There is no phase separation, i.e, discontinuities in n{^), 
for J/t — > since the ji—n relation is that of the underlying 
tight binding system and free of any singularity. The phase 
transitions with changing n are all second order. With 
growing J/t, however, some phase boundaries become first 
order and regimes of PS will emerge. 

(ii) Strong coupling: For J/t — > oo, it makes sense to 
quantise the fermion spin at site in the direction of 
the core spin Si, and project out the 'high energy' un- 
favourable state. This leads to an effective spinless fermion 
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Fig. 3: Colour online: The magnetic ground state in the RKKY 
limit, showing the dependence of the ordering wavevector Q 
on electron density and particle-hole asymmetry (via t'). The 
legend for the various states is shown on the right. The calcu- 
lations were done at weak coupling, J = 0.5. Note the growing 
asymmetry of the phases (about n = 1) as t' increases. It is also 
clear that if the hopping parameter t' changes (due to pressure, 
etc) the magnetic ground state can change even if the electron 
density remains fixed, as discussed for 4/ systems in [20] . This 
is particularly prominent in the top right hand corner of the 
figure. In constructing this phase diagram we have ignored a 
narrow sliver of phase separation near n — 1. 

problem whose bandwidth is controlled by the average spin 
overlap (Si.Sj) between neighbouring sites. The overlap is 
largest for a fully polarised state, and the FM turns out to 
be the ground state at all n 7^ 1. At n = 1 'real hopping' 
is forbidden so the fermions prefer a G type AFM back- 
ground to gain kinetic energy 0{t'^ / J) via virtual hops. 

The FM and G type AFM have a first order transition 
between them with a window of phase separation, eas- 
ily estimated at large J/t. The fully polarised FM phase 
has a density of states (DOS) which is simply two 3D 
tight binding DOS with splitting J between the band cen- 
ters. If we denote this DOS as NpMii^, J) then the en- 
ergy of the FM phase is SpAiif^, J) = Jt^ac J)Luduj, 
and the particle density is n{ii,J) = ^^^Nfm{^i J)dijJ. 
There will be corresponding expressions when we con- 
sider electrons in the {tt, tt, tt} AFM background, with 
DOS Nafm{^^, J)- Once we know ^ = i^AFM ^^'"^^ 
isfies £pi\i(fi,J) = £AFM{^^^J) we can determine the 
PS window from the density equations. Since the FM 
phase has a dispersion e^^^ = eo.k ± J I'^i where eo.k = 
—2t{coskxa + coskyU + coskza), and the AFM phase has 

dispersion e^^^^ = ± y^e^^~f~(J/2)^, it is elementary to 

work out IJ-AFAI- The analysis can be extended to several 
competing phases. It is significant that even at J/t ~ 10, 
which might occur for strong Hunds' coupling in some ma- 
terials, the FM phase occurs only between n = {0, 0.7}. 

(iii) Intermediate coupling: The intermediate coupling 
regime is where one is outside the RKKY window, but not 
so large a coupling that only the FM and G type AFM are 
possible. Towards the weak coupling end it implies that 



the planar spirals begin to pick up a net magnetisation, a, 
and now become 'conical' phases. Windows of phase sep- 
aration also appear, particularly prominent between the 
{Qx, Ti", tt} and G type AFM (near n = 1), and suggest the 
possibility of inhomogeneous states, etc, in the presence of 
disorder. The prime signature of 'physics beyond RKKY', 
however, is that the RKKY planar spirals now pick up a 
net magnetisation and much of the phase diagram starts 
to evolve towards the ferromagnetic state. 

(iv) Finite temperature: The TCA based MC readily 
captures the FM and {7r,7r,7r} AFM phases at all cou- 
pling. However, it has difficulty in capturing the more 
complex spiral. A, and C type phases when we 'cool' from 
the paramagnetic phase. In the intermediate J regime it 
usually yields a 'glassy' phase with the structure factor 
having weight distributed over all q. In our undertstand- 
ing this is a limitation of the small cluster based TCA, 
and the energies yielded by VC are better than that of 
'unordered' states obtained via MC. To get a feel for the 
ordering temperature we have calculated the energy differ- 
ence A£, defined earlier, as often done in electronic struc- 
ture calculations. This provides the trend in Tc across 
the phases. Fig. 2, and wherever possible we have included 
data about actual Tc (symbols) obtained from the MC cal- 
culation. Broadly, with increasing J the A£ and Tc scales 
increase but the number of phases decrease. The Tc of the 
G type AFM is expected to fall at large J but even at 
J/t ^10 it is larger than the peak FM Tc. 

(v) Interplay of FS and coupling effects: Till now we 
have looked at the particle-hole symmetric case where 
t' = 0. The tight-binding parametrisation of the ab ini- 
tio electronic structure of any material usually requires a 
finite t', in addition, possibly, to multiple bands. We will 
use the t — t' parametrisation of band structure due to its 
simplicity. It will also allow us to mimic the physics in the 
4/ metals. 

At weak coupling the magnetic order is controlled as 
usual by the band susceptibility, xo(Qj ti) which, now, also 
depends on t' . At fixed n the magnetic order can change 
simply due to changes in the underlying electronic struc- 
ture. Our Fig. 3 illustrates this dependence, where we use 
J ~ 0.5 to stay in the RKKY regime and explore the vari- 
ation of magnetic order with n and t' . The range of t' 
variation is modest, ~ {0 — 0.3}, but can lead to phase 
changes (at fixed n) in some density windows. We have 
cross checked the phases with the peak in Xo(q)- 

A complicated and more realistic version of this has 
been demonstrated recently [20] in the 4/ family for the 
heavy rare earths from Gd to Tm. These elements all 
have the same hep crystal structure, and the same conduc- 
tion electron count, 5d^6s^, so nominally the same band 
filling. However, the electronic structure and Fermi sur- 
face changes due to variation in the lattice parameters 
and unit cell volume (lanthanide contraction) across the 
series. It has been argued [20] that this changes the loca- 
tion (imax of the peak in xoi a-nd explains the change in 
magnetic order from planar spiral (in Tm) to ferromag- 
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Fig. 4: Colour online: Top: the magnetic ground state at 
n = 1.7 for varying t' and J. The ordering wavevector is 
marked on the phases. The magnetic order has a pronounced 
dependence on both the 'bandstructure' (through t') and the 
electron-spin coupling. We highlight three kinds of parame- 
ter variation, (i) Varying J at fixed t' , the points on the y 
axis, shows how changing electron-spin coupling can change 
the ground state, (ii) Varying t' at weak coupling, J = 1, il- 
lustrates how bandstructure affects the RKKY magnetic order, 
(iii) In the 4/ elements we think what happens is a combina- 
tion of (i) and (ii) above, as shown by points on the diagonal. 
Bottom: An impression of the real space spin configuration 
for the three parameter sets (i)-(iii) in the top panel. Each 2x2 
pattern is for at', J combination. The bottom left spin in each 
pattern is set on the reference site R = {0, 0}, say. The neigh- 
bouring three spins are at {a;,0}, {0,j/}, and at {x,y}, where 
X and y are unit vectors on the lattice. There is no variation 
in the z direction so we only show the in-plane pattern. Top 
row: scan (i) above, changing J at t' — 0. Middle row: scan 
(ii), changing t' at J = 1. Bottom row: scan (iii), simultaneous 
change in t' and J. 



netism in Gd. A similar effect is visible in our Fig. 3 where 
at n = 1.7, say, the ordering wavevector changes from a 
spiral to FM as t' changes from zero to 0.15. In this sce- 
nario, J does not affect the magnetic order but merely sets 
the scale for Tc- The RKKY interaction strength scales as 
J^ff ~ J'^S{S+1), and a similar scaling of the experimen- 



tally measured Tc is taken as 'confirmation' of the RKKY 
picture. 

Should'nt we also worry about the effect of the grow- 
ing Jeff(S) on the magnetic order itself! If the maximum 
Jeff, for Gd with 5* = 7/2, were smaller than the effective 
hopping scale t, then wc need not - the RKKY scheme 
would be valid for the entire 4/ family. However, mea- 
surements and electronic structure calculations [3] in Gd 
suggest that J - 0.3eV and Jeff {7/2) - IcV. The effec- 
tive t is more ambiguous, since there are multiple bands 
crossing the Fermi level, but the typical value is ~ 0.3eV. 
This suggests Jeff/i ^ 3, clearly outside the RKKY win- 
dow! What is the consequence for magnetic order, and 
physical properties as a whole? 

Fig. 4 shows the t' — J magnetic phase diagram at T = 
for n = 1.7. At t' = 0, the vertical scam, changing J reveals 
how the ordered state changes with increasing J even with 
electronic parameters (and hence xo and FS) fixed. We 
have already seen this in Fig. 1 The spirit of RKKY is to 
assume J ^ 0, and move horizontally, changing t' across 
the series so that one evolves from a planar spiral to a 
ferromagnet. We suggest that in the / metals, the param- 
eter points are actually on a 'diagonal', with increasing 
t' (our version of changing electronic structure) being ac- 
companied by increase in Jeff- To capture the trend we 
set, t' = and Jeff = 1-0 for S* = 1, where the system 
is known to be a spiral, and i' = 0.1 and Jeff = 3.0 for 
S = 7 /2 (the case of Gd), and explore the linear varia- 
tion shown in Fig. 4. This paramctrisation is only meant 
to highlight the qualitative effect of changing electronic 
structure and Je// and since real t' values, etc, would need 
to be calculated from an ah initio solution. 

Within this framework, while the small S result is same 
for both RKKY and explicit inclusion of Jeff, the order 
obtained at intermediate S depends on whether one ig- 
nores Jeff (as in RKKY) or retains its effect. For a given 
t' the phase on the diagonal is quite different from the 
phase on the horizontal line. 

In fact there is evidence from earlier ah initio calcula- 
tions [21] that in addition to unit cell volume and c/a ratio, 
the strength of the 4/ moment (and so Jeff) also affects 
the magnetic order. As an illustrative case, the optimal 
spiral wavevector in Ho evolves towards Q = {0,0,0} as 
the effective moment is (artificially) varied from 2^b to 
AfjLB (Fig. 2 in Nordstrom and Mavromaras [21]). If mag- 
netism in this element, and the 4/ family in general, were 
completely determined by RKKY there would be no de- 
pendence on Jef f ■ In fact the authors suggested that one 
should re-examine the basic assumptions of the 'standard 
model' of 4/ magnetism [5], which gives primacy to the 
RKKY interaction (and magnetoelastic effects) since the 
ah initio results suggest a role for the effective exchange in 
the magnetic order. Our aim here has been to clarify the 
physics underlying such an effect within a minimal model 
Hamiltonian. This approach would be useful to handle 
non coUinear phases in complex many band systems, with- 
out any weak coupling assumption, once a tight binding 
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parametrisation of the electronic structure is available 

Let us conclude. We have examined the Kondo lattice 
model with large S spins and established the ground state 
all the way from the RKKY regime to the strong coupling 
limit. The intermediate coupling window reveals a com- 
petition between RKKY effects, which tend to generate a 
planar spiral, and the tendency to gain exchange energy 
via ferromagnetic polarisation. This generally leads to a 
'conical' helix, giving way at strong coupling to the double 
exchange ferromagnet. Using these results we rc- visited 
the classic 4/ magnets to demonstrate how the magnetic 
phases there arc probably controlled non RKKY spin- 
fermion effects. One can add anisotropics and magneto- 
elastic couplings to our model to construct a more com- 
prehensive description of 4/ magnetism. 

We acknowledge use of the Beowulf cluster at HRI, and 
thank Sanjeev Kumar and B. P. Sckhar for collaboration 
on an earlier version of this problem. 
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